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Abstract. We investigate the behaviour of a chain of interacting Brownian par- 
f— ^ ', tides with one end fixed and the other moving away at slow speed e > 0, in the 

CNJ ' hmit of small noise. The interaction between particles is through a pairwise poten- 

tial U with finite range 6 > 0. We consider both overdamped and underdamped 
dynamics. 
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'^ I The behaviour of a Brownian particle moving in a potential well and acted upon by a 

^2 ' linearly increasing force is widely used to model the mechanical failure of molecular bonds 

arising in dynamic force spectroscopy experiments |LCST07l IDFKU03| lREB"'"06l IFriOSj . 
This began with the work of Bell |Bel78] and was developed further by Evans and Ritchie 
[ER97]. 

K*" I Let Qs denote the length at time s of a bond that is fixed at one end and has a harmonic 

^^ ■ spring attached to the other. If the spring moves linearly at speed e > 0, the motion of Qs 

f^ , is typically modelled according to an SDE of the form 

^ ■ dqs = i-U'iqs) + es) ds + a dWs , 

Q ■ where U{q) denotes the bond energy (e.g. Lennard- Jones potential), Wg is a standard 

^D . Brownian motion and o" > is the (small) noise intensity. Note that this model assumes 

the motion is overdamped. Rupture of the bond corresponds to the first time Qs escapes 
from the stable well of the effective, time-dependent potential, H{q,es) = U{q) — esq. 
^ , The effect of the external force is to lower the barrier height of H, thus making escape 

^ '_ more likely. The main objective is to study the distribution of first-breaking times and 

how the mean first-breaking time scales with the pulling speed e. Typically, two pulling 
speed regimes are considered. 

For very slow pulling, the particle is able to escape the well through a large deviation 
event before the potential has changed significantly and the energy barrier is still large. In 
order to apply the standard theory valid for time-independent potentials [Kra401 Eyr35 



C^ 



IBEGK041 IFW98| . the adiabatic approximation is used: at any given time s, the bond 
has an instantaneous rate of rupture, k{s), and the probability of survival until time s, 
denoted P{s), decays according to P{s) = —k{s)P{s). Note that l/k{s) is the usual 
Eyring-Kramers formula |Kra401 Eyr35[ IBEGK04] applied to H at time s. 



As the speed of pulling increases, the energy barrier at the time of rupture becomes 
smaller. If pulling is sufficiently fast, the barrier may be close to vanishing completely 
when rupture occurs. This means that the external force, given by es, is almost equal to 
the maximum slope of U, which occurs at the point of inflection between its minimum and 
maximum, i.e. maximum slope is U'{cq), where U"{cq) = 0. For times s at which es is 
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close to this critical force, the effective potential H is almost cubic near its minimum. This 
leads to a different rupture rate than that given above, although still calculated within 
the Kramers framework. 

It is interesting to consider what happens as the pulling speed increases yet further 
and the Eyring-Kramers formula is no longer applicable, nor the adiabatic approximation 
underpinning the above approach. In this paper, we consider this situation in a model 
related to that above. More precisely, we consider a chain of two identical bonds in series 
with one end fixed and the other being pulled at a constant rate e. Both overdamped and 
underdamped dynamics are treated. We are interested in which of the two bonds breaks 
first and how this depends on e and the noise intensity a. As above, the dynamics near 
the inflection point of the bond energy U play an important role and will be the focus of 
our analysis. Roughly, we flnd that for e > a^'^, the right-hand bond breaks first, while 
for e < o"^'^, both have an equal probability of breaking in the limit of small noise. Thus 
e = o"^'^ represents the threshold at which the adiabatic approximation becomes valid. 

To our best knowledge, the first work to tackle rigorously such models of bonds under 
an external, time-dependent force was [AB09| . There the authors consider a similar model 
of two bonds in series as above, but with an additional assumption that U is cut-off strictly 
convex. The breaking event corresponds to the first time one of the two bonds exceeds the 
range of U. Roughly speaking, it is shown that for e > a, the chain always breaks on the 
right-hand side, whereas for e < a, each bond has an equal chance to break in the small 
noise limit. Thus the threshold between the different types of behaviour is different from 
that found in the present work, where the bond energy U is taken to be smooth (but also 
with finite range). In principle, the results of |AB09| can be extended to arbitrarily many 
bonds in series j All 10 . 



The behaviour of several bonds in series has also been considered by many authors, for 
both time-dependent and -independent external forces. The situation when the external 
force is constant, i.e. one initially stretches the chain by some amount and then fixes 
both endpoints, has been considered for harmonic potentials [Lee09] and Lennard-Jones 
potentials [SDGOGj . In the harmonic case, it is shown analytically and numerically that 
the probability to break at either endpoint is half that of breaking at any non-extremal 
point, which all have the same probability. In the Lennard-Jones case, the motion is not 
assumed to be overdamped, i.e. the authors consider the equation 

r)J-f 

mqi{s) = -mjqiis) - 7^(g(s)) + ^/2^W^(.s) , 2 < i ^ iV - 1 , 

where the Wi are independent Brownian motions, q{s) = (9i(s))i=i G ^ is the chain 
configuration, with qi and q^ constant, and H is the total potential energy of the chain. 
It is assumed beforehand that one bond is close to breaking and all others are close to 
minimal energy so that a quadratic approximation can be used for H. Then the breakage 
rate is calculated, alongside simulations, using a multi-dimensional version of Kramers' 
theory developed by Langer |Lan69j . These show that the breakage rate is lower the 
closer the chosen weak bond is to the chain endpoints. However, the simulations also 
show that the harmonic approximation for H may fail for bonds near the breaking bond, 
which the authors there suggest may explain some discrepancies between the theory and 
simulations. As they point out, Langer's theory, as well as the classical Kramers theory 
for a single particle, requires a harmonic approximation for H. 

The case of a chain with one end fixed and a linearly increasing force applied at the 
other end has been considered by Fugmann and Sokolov in [ FS09b. TS09a] to model the 
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mechanical failure of a polymer chain. More precisely, they consider the vector q(s) = 
{(li{^))iLo ^ IK^"''"'^, with go = 0, which evolves according to the SDE 

f)J-f 

dqi{s) = -j—{q{s),es) ds + a dWi{s) , 1 ^ i ^ TV , 

where H{q,es) is the time-dependent potential energy of the chain, given by 

H{q, es) = J2 ^(9« - 1j) -^sqN , 
^ i<j ^ N 

and U is the Morse potential. A break is said to occur when q{s) overcomes an energy 
barrier of the effective potential, H. Numerically, they show that for a high pulling speed, 
e, only the right half of the chain contributes to the breaking event and the probability 
increases as you move towards the right. For smaller e, they show that the breakpoint is 
more uniformly distributed along the whole chain. In their analysis, they assume that the 
rupture dynamics of different bonds are independent and then apply the one-dimensional 
theory. The only thing left to do then is to analyse how much force each bond feels, which 
depends on the pulling speed. For lower pulling speeds, it is assumed that each bond feels 
the same force, whereas for higher pulling speeds the chain is approximated by a harmonic 
chain. 

This paper is organised as follows. In Section O we introduce our model and deduce 
equation (j2.2p . which is our main object of study. Our results are then stated in Theorems 
12. H 12.21 and 12.31 In Sections El [4] and [5] we give the proofs. 

The following notation is used in this paper: 

• By xi{t,e) X X2{t,e) we mean that for two functions xi(t,e), X2{t,e), defined for 
t in an interval I and < e ^ Eq, there exist constants c± > such that 

c-X2{t,e) ^ xi{t,e) ^ c+X2{t,e) 

for all t e / and < e ^ eo- 

• By Ox{l) we mean that lim^_!.o 02^(1) = 0. 

• xi < X2 means that there exists a constant c > such that cxi ^ X2 for xi, X2 > 
sufficiently small. 

• xi <C X2 means that xi = o{x2). 

• We shall write f*-0'io to denote probability conditioned on the relevant process 
starting at time to in position go- 
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2. The model and main results 

Three particles, qL, q and q^ in M, interact with each other via a pairwise potential U. 
We assume that U is smooth with finite range 6 > and a unique minimum at < a < 6, 
with U"{a) > 0. We also assume that there is a unique cq G {cL,b) such that U"{co) = 0. 
The particle qi is fixed at the origin and the position of g/j at time s ^ is given by 
Qr{s) = 2a(l + £s), where e > is a small parameter. We study the behaviour of the 
middle particle, with position at time s given by qg. Initially, it has position qo = a so that 
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the distance between neighbouring particles is a. The middle particle evolves according 
to an SDE of the form 

dqs =psds , 

BJ-f 

e^'^dps = -psds- —{ps,qsj£s)ds + adWs , 

where Ws is a standard one-dimensional Brownian motion with Wq = 0, cj > is the noise 
intensity, /? E M and H{p,q,es) is given by 

2 

H{p, q, es) = Y + U{q) + C/(2a(l + es) - q) . 
Rescaling time as t = es, this is the same in law as solving 



1 , 

aqt = -Pt dt 

e 



1 . IdH 



a 



e^~'dpt = — Pt dt — (pt, qt, t)dt + ^ dWt 



(2.1) 



t ■ 



e e oq \Je 

The length of the chain {qL,Q,QR) increases linearly with t. Clearly, if we wait a long 
enough time, the distance between q and at least one of its neighbours must become 
greater than the range of U. In this case, these particles no longer interact and the chain 
can be considered broken. Since U has a minimum at o, it is energetically preferable for q 
to move towards either qi or q^. Letting e = e{a), our aim is to determine how the speed 
of pulling affects which of these two possibilities occurs in the limit as o" ^ 0. 

We easily check that the configuration of equally spaced particles satisfies dgH = 
0, dqH > until time to, where a(l + to) = cq. Thus until this time it is a stable 
configuration and so we expect qt^ ~ a(l + to)- For t > to, this configuration becomes 
unstable and new minima emerge. So we expect qt to quickly move away from the chain 
midpoint and towards one of these newly formed minima. Note that as a function of q, 
H is symmetric about q = a(l + 1), but its time-dependence introduces asymmetry as we 
shall see. Once qt has approached one of these new minima, we expect it to stay there 
as the energy barrier to escape becomes higher. The evolution of the chain, therefore, is 
determined by its behaviour around the bifurcation of -ff at t = to, which we shall now 
consider. 

Letting zt = a(l + 1) — qt^ we express the term dH/dq appearing in (12. ID in terms of z. 
By a Taylor expansion in space, we find 

BM 

-Q^iPuqt,t) = U'iqt) - U'{2a{l + t) - qt) 

= U'{a{l + t) - Zt) - U'ia{l +t) + Zt) 
^ -2U"{a{l + t))zt - ^U^^Hail + t))zf . 

Assuming that there is < T < to such that for t G [to — T, to+T], U^^\a{l+t)) is negative 
and bounded away from zero (see comment below), we have by a Taylor expansion in time, 

-2;7"(a(l + t))zt - -t/(^)(a(l + t))zf ^ 2a{t - to)zt - Czf . 

We remark that this assumption about [/'•^^(a(l + t)) should not have much effect. At 
most, the right-hand side above would have a +Czf term appearing, but in either case 
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this term is very small for z and t — to close to zero, which is where most of our analysis 
will take place. 

Making the space and time transformations q = a{l + 1) — q, p = a — p/e and t = t — tQ, 
as well as normalising constants to one, we arrive at the SDE 

dqt =Ptdt , 

ef^dpt = -pt dt + -{tqt - qf + e) dt + ^ dWt . ^^'^^ 

This will be our main equation for the rest of this paper. By the above discussion, 
understanding how its solution behaves will be a good indication of the behaviour of 
the original chain. Equation (|2.2p represents the motion of the particle q in the potential 
(l/e)y((7, t) := {l/e){—^tq'^ + jq'^) with an additional +1 force giving the particle a small 
bias towards the right. This force comes from pulling the chain (g^, q, g/j) and corresponds 
to the fact that in the absence of noise, q does not just stay at the chain midpoint a{l + t), 
but lags behind by a small amount. 

Rephrasing the discussion after ()2.ip . the function V represents the energy of a given 
chain configuration. For negative times, the origin, corresponding to equally spaced par- 
ticles, minimises V. When t = 0, V undergoes a symmetric pitchfork bifurcation at the 
origin. For positive times, V has two minima located at ib\/t. For t > large enough 
these minima at ib-v/t correspond to the configurations where q is a distance a from q^ or 
qn, respectively, and more than b from the other. In terms of (j2.2p . the aim of this paper 
can be roughly stated as to determine whether q moves towards +^/i oi —^/i as t becomes 
positive, which corresponds to the chain 'breaking', and how it is affected by the speed of 
pulling. 

There are several ways in which one may rigorously define the chain to break. One 
possible definition is that the chain breaks as soon as the distance between q and one 
of its neighbours exceeds the range of the pairwise potential U. Then the chain either 
breaks on the right- or left-hand side, depending on whether qji — q > b 01 q — qL > b, 
respectively. This was used in |AB09] . Alternatively, one may consider the chain to break 
as soon as the chain configuration reaches a neighbourhood of one of the energy minima 
that emerges after the bifurcation. In the above formalism, this means the process qt 
reaching a neighbourhood of ib\/t. We shall avoid making this choice by instead giving 
a precise description of the behaviour of qt that contains more information than any of 
these possible definitions. 

Equation (j2.2p in full is not something we can treat. But there are two obvious simpli- 
fications: the first is to omit the qf term in the equation for p, leading to a linear equation 
that can be solved explicitly; the second is to neglect mass, taking e"dp = 0, and consider 
the overdamped equation. We treat both of these and obtain satisfactory results. 

Taking a = e"+i/2 for a > —1/2, we firstly consider the linear SDE 

d(??=P?dt, 
e^dp^t = -p° dt + - {tq^ + e)dt + e" dWt . 

Denoting by P'' the law of the solution with vanishing initial condition at time s < 0, we 
have the following result: 

Theorem 2.1. Let q^ be the solution of h2.3\) . If a > 1/4 then 

lim lim inf P'* < lim q^ = +c« > = 1 , 

elO s—>-—oo It—^oo J 
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while if a < 1/4 then 

limliminf P"* | lim q'^ = +00 j = lim lim inf P'' ( lim q^ = -00 j = 1/2 . 

£4-0 s—>-—co li— >oo J eIO s— !>— 00 It— s>oo J 

This theorem shows that the threshold between fast and slow pulling regimes is given 
by a = 1/4 and is independent of (3. However, we also note that if we start the processes 
at a finite negative time — T, then neglecting mass does have an effect, but only for /3 < 0. 
Indeed, for —1 < /3 < and zero initial conditions, the threshold becomes a = (1 + /3)/4 
(see [All 10] ). This result has some clear limitations: for t > 0, q^ shoots off quickly to 
±00 as the drift becomes ever more repelling, while the solution of (j2.2p is prevented from 
doing this by the nonlinear term and so is more likely to return to the origin. 

Secondly, we neglect the mass term e^dp in ()2.2p and consider the one-dimensional 
overdamped equation 

dqt = - {tqt - qf + e) dt + ^ dWt . (2.4) 

Again letting P* denotes the law of the solution with vanishing initial condition at time 
s < 0, we have 

Theorem 2.2. Let qt solve \2.4^ . There exist constants ci,7 > such that if ti = 
ci-^e| Ino"! then 

(1) (Fast Pulling) for any cj^/^| Ino-p/^ < £{a) <C 1, 

limliminf P'' i inf -^ > 7 I = 1 . 



(2) (Slow Pulling) for any cr^l Incrp < e{a) < a^/'^\ hia\~'^'^l^, 





and 



lim lim sup 

cJ-O s-s>-oo 



lim lim sup 



mt ^ > 7 
ti^t ^ 



qt ^ \ I 

t?lvi^~'i~2 



0. 



Letting a = £"^+1/2 above gives a = 1/4 as the threshold between the different regimes, 
as found in Theorem 12.11 We also remark that the threshold between fast and slow pulling 
regimes here differs from that in [AB09| . where it is roughly e = a. This difference can be 
attributed to the bifurcation of V: for t near zero, V is almost flat and so the additional 
+1 force requires slower pulling, or stronger noise, to be counteracted than it does in 
|AB09] ■ where the potential has positive curvature bounded away from zero. 

Having considered two simplifications of (|2.2p , we finally return to the full solution itself. 
Intuitively, by taking /3 large, the effect of the mass term e^ should become small and the 
solution should behave like that of the overdamped equation ()2.4p . So if we show that for 
suitably large /3, the difference between the two solution stays small, we can use Theorem 
to tell us about (j2.2p . This leads us to the following result: 



Theorem 2.3. Let qt solve h2. 2\) with (3 > 2. There exist constants ci,7 > 0, independent 
of j3 and a, such that for ti = Ci-y/e| Ino"! and any t2 > ti, 

(1) (Fast Pulling) if a'^^'^l Inap/^ < £{a) < 1 then 
lim lim inf P'* < inf —= > 7 

0-40 s-S>-oo [ti ^ i ^ i2 ^/t 
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(2) (Slow Pulling) while if < 6 < (3/2-1 and cj2/(i+25) ^ ^^^^ ^ ^4/3 1 ij^^|- 13/6 
then 



lim lim sup 

o"4-0 s_^_oo 



inf ^>-/}- 1/2 







and 

1 «4 I 



lim lim sup 

o"4-0 s^— 00 



sup ^ < -7 >> - 1/2 

ii < t < t2 VI 



Note that we only consider finite time intervals here and that there is a slight difference 
between the lower bound on e in (2) above and in Theorem 12.2^ 2). This second point 
is related to the fact that the mass is of the form e^ . However, it does not affect the 
threshold between the two regimes. 

We finally note that (j2.2p is not suitable for considering the chain 'breaking' due to a 
large deviation event. In that case, our expansion of the potential is not valid. 

3. The linear model 

In this section, we give the proof of Theorem 12 . 1 1 and are therefore dealing with equation 
(j2.3p . This equation is simple enough to have an explicit solution in terms of Airy functions 
Ai(t), Bi(t) (see |dlm07j ) : using the fact that both Ai{z) and Bi(z) solve the equation 
w"{z) = zw{z), and that the Wronskian Ai'(t)Bi(t) — Bi'(t)Ai(t) = l/vr, it can be checked 
that the process 

g°(t) =7re(i-2^)/3 f-Ai{t{e, /3)) / q-^^'-''^''^ Bi(s(e, /3))(ds + e^dWs) 



+ Bi(i(e, /?)) / e-^^*-^)^'" Ai(s(e, /3))(ds + e^dVF,)) 
J —00 / 



(3.1) 



is the (almost surely unique) solution of (|2.3p with zero initial conditions and starting 
time sent to —00. Above, we have set s(e,/?) = e~^^~^^''^{s + e^~'^/4). The semi-infinite 
stochastic integrals are to be understood as limits of finite stochastic integrals from —T to 
t as T — )• 00. The existence of a limit process q^{t) of this procedure, for all t in intervals 
[— 5, 00), S" > 0, is easy to see in the present case: Clearly, when —00 is replaced by 
— T, (|3.ip defines a Gaussian Markov process on [—T, 00), with mean and variance readily 
calculated (albeit with long formulas, which is why we do not display them here). By 
using the asymptotic expansions 

Ai(s) = -^s-'/^ e-2-''V3 (1 + 0(5-2/3)) ^ 5 > 1 

2^/^T 



Bi(s) = ^s- 1/4 e2^^/V3 (1 + o(s-2/3)) , s > 1 

Ai(s) = ^|s|-i/^cos(2|s|3/2/3 - 7r/4)(l + 0{\s\'^/^)) 
Bi(s) = ^\s\-^/^ sin(2|sp/2/3 - ^/4)(1 + 0{\s\-^/^)) 





(3.2) 




(3.3) 


s< -1 


(3.4) 


s< -1 


(3.5) 



it is clear that the distribution of qt at any time —S converges to a non-degenerate Gaussian 
as T — )• C!0, and thus (13.11) is well-defined. 
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The behaviour of g*^ as t — >• oo can be described in a straightforward way by considering 
the 'renormahzed process' 



^W = ^£(i-2/3)/3Bi(t(e,/3))'^ ^^> 



Indeed, we then have 



Proposition 3.1. lirat^oo Qt exists almost surely, and is a Gaussian random variable with 
mean 

m = e(i+/5)/3e-i^^""' (3.6) 

and variance 

V = e2"+(i+/3)/3 e-i--^^ r e-'^-^^'^^ M{sf ds . (3.7) 

J — QO 

Proof. Let us first investigate the deterministic integrals in ()3.ip . Since Bi(s) < Bi(t) for 
s < t and t ^ 0, we have that the deterministic integral in the first line of ()3.ip . after 
renormalisation, is bounded by Ai(t(e,/3)) J_^ e**^ '^ ds for large enough t, and thus 

converges to zero as t — )• oo due to ()3.2p . It is known [dlmOTj that J^ e^*Ai(s) ds = e^ '^ 
for all p > 0, and thus the limit of the corresponding (renormahzed) integral in the second 
line of (|3.ip is given by 

oo /'OO 



et^ ^ Ai(s(e, /?)) ds = e^^+'^^Z^ e-p ' / e^^^" ^" Ai(s) ds = e^'+i')^-" ^-^2 

oo J —oo 

Clearly, this is also the limit of E[q{t)) as t — )• oo. 
The stochastic integrals in q{t) are given by 

Ji(t) = e^h{t) f Ms) dWs , J2{t) = 8" f h{s) dWs , 

J —oo J —oo 

with 

h{t) = ^r^rik^ /i(-)=Bi(.(£,/3))e--^/2, /2(s) = Ai(s(£,/3))e--'/2. 
Bi(t(e,/?)) 

By the time change Ti(f) = e^°' f_^ /i (■?) ds, Ji(t) equals h{t)B^^(j.-^ in distribution, where 
Bs is a standard Brownian motion. By the law of the iterated logarithm, 

hm sup , ''^'^^^' = 1 (3.8) 

t^oo h{t)yj2Ti{t)\n\nTi{t) 

almost surely. Again using Bi(s) < Bi(f) for s <t and t ^ 0, we find 

rt . \ 1/2 



h{t)^/M^ ^ Ai(t(e, /3))e° ( / e"^" ds 

\J —oo 



which by (j3.2p converges to zero superexponentially fast. By (13. 3p it is easy to see that 
Inlnri(t) grows only proportionally to Int, and thus the denominator on the left-hand 
side of (j3.8p converges to zero. It follows that Ji(t) — )• as t — )• oo almost surely. J2, on 
the other hand, is a square-integrable martingale, and thus converges almost surely. Each 
J2{t) is Gaussian, and thus so is the limit. It has mean zero and variance 

/oo /'OO 

fl{s) ds = e^" / e^""'^ Ai(s(e, I3)f ds . 
-00 J —00 

The same change of variable that was employed to get (|3.6p yields (j3.7p . D 
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Since Bi(t)e~* diverges as t — >• oo, Proposition 13.11 means lim(_!.oo k'^COl = co almost 
surely. Whether the divergence is to plus or minus infinity is determined by the sign of 
goo = linit-5>oo Qt- Qoc is a Gaussian random variable with mean m = m(e) > given 
by (|3.6p . and variance v = v{e) given by (|3.7p . Thus if lime_^o "i(e)/\/'^(£) = oo, then 
lim(;_j.o lP(hmt^oo Q^it) = +oo) = 1, as the distribution of goo concentrates on the positive 
half line. On the other hand, if lim^-^o "i(£)/\/^(^) = 0) then lime^oP(li™t^oo <?"(*) = 
+oo) = 1/2, as the distribution of goo becomes spread out and P(goo > 0) ^ 1/2 as e ^ 0. 

We now determine the circumstances under which each of the above cases occurs. Define 



/oo 
-oo 



We have 



Lemma 3.2. There exist constants ci and C2 such that 
(i) limp^ooP^/'e-2^V3j(p) = ci, 
(ii) limp^o P^^'^ e~^P /^ J{p) = ci . 

Proof. Consider first the case p — t- oo. Then, 

e^P' Ai2 (s) ds pi/2 e- V/3 ^ ^ ^-2pV3+2p ^1/2 



as p ^ 00. For 5 > 1 we use (|3.2p and find 

e^P' Ai\s) ds = ^ g-4«3/2/3+2ps ^-1/2(1 ^ 0(s-3/2)) ds 

00 

^-p3(4t3/2/3-2t) ^-l/2(-^ ^ 0(p-3t-3/2)) dt , 

l/p2 

where we used the substitution s = p^t. Decompose the integral as j'y 2 = /i/ 2 + /i/4- 

The first of these is bounded by Ce^ '^ for some C > and we can ignore it. For the 
second, we have 0{p~^t~^i'^) = 0{p~^) and can take this outside the integral. Then by 
the Laplace method, 

00 , /'OO 

p3(4t3/2/3-2t) ^-1/2 ^^ ^ g2p3/3 / g-p3[(i_i)2+o(j_i)3] ^_i/2 ^^ 



e 

1/4 ^1/4 

Thus (i) holds with ci = j^. For (ii), we use that 

/oo /"OO 

e2P"Ai2(s)ds^ / Ai2(s)ds = const 

as p — > 0. With ()3.4p we then get 

/ e^P'A\^{s)ds = - f e^P'\s\-^/'^cos\2\sf/^/3-Tr/A){l + 0{\s\-^/^))di 

J ~00 " J—QO 
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where in the last hne we used the substitution t = —ps. As p — t- 0, the integral in the last 
line above converges to 

2 70 2^/2 ' 

which proves (ii) with C2 = ^Ly^ . D 

When substituting p = e(i~2/3)/3^2 into the last lemma, we find that 

Thus, m^e)/ yjv{e) = C'e"""*'^'^, which completes the proof of Theorem 12.11 

4. The Overdamped Model 

The proof of Theorem 12.21 consists of two parts. First, we show that given an arbitrary 
compact set X containing a neighbourhood of the origin (say X = [—1,1]), and given 
an arbitrary negative time — T, the process qt starting at at time — cxd will be in X 
with very high probability at time —T. In a second part, we then take advantage of the 
Markovian nature of the process to restart it at time — T and to show that the conclusion 
of Theorem 12.21 holds uniformly over initial conditions belonging to X at time —T. These 
two steps are formulated as Proposition 14.11 and Proposition 14.21 below. Theorem 12.21 is 
then an immediate consequence of these two results. 

The first part can be formulated as follows: 

Proposition 4.1. Let X = [—1, 1], let T ^ 1, and let qt solve ^2.4\ )- Then, we have 

lim liminfP^(g_T £ X) = I . 

o-,e— >0 s— >— oo 

Proof. Applying Ito's formula to the function q^^ q^ , we obtain 

-E(?2 = E -g2 _ -qf + 2gj + — U - -Eg2 + e + — , 

at \ £ £ £ J £ £ 

where we made use of the fact that t ^ — 1. The claim then follows at once from the fact 
that we assumed that e — t- and cr — t- 0. D 

The remainder of this section is devoted to the proof of the following statement, where 
we denote by P~^'^ the law of (12. 4p with initial condition q^x = x € X. 

Proposition 4.2. Let qt solve {2.4\)- There exist constants ci,7 > such that if ti = 



c\\l £\ Ino"! then 

(1) (Fast Pulling) for any a^/^\ lno-p/3 < s{a) <C 1, 

hm inf P-^'^ I inf i^ > 7 I = 1 . 

cr^OxeA' [tlsit^/t J 

(2) (Slow Pulling) for any cr^j Incrl^ < £{a) < cr^/^l lna\^^^/^, 





lim sup 

0-40 x&x 






and 

r^ \ ri-i 1 



lim sup 



-r,x J sup ^ < _^ [, _ 1/2 
M ^ t yt 
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Our approach in this section is based on that developed by Berglund and Gentz in 
[BG021IBG06J . They consider similar equations to (j2.4p . but with drift terms {l/e)f{q,t) 
such that f{q,t) = —f{—q,t) and /(0, 0) = 5^/(0,0) = 0. A simple example of such an / 
is f{q,t) = tq — q'^. Our additional drift term arising from pulling means that we cannot 
directly apply their results, except in a few cases as will be made clear. 

4.1. Fast pulling. We begin by considering the fast pulling regime from Proposition 14.21 
In this case, the noise in the systein is not strong enough to overcome the asymmetry 
caused by pulling and the qualitative behaviour of qt is the same as that of the deterministic 
solution qf""^ of the ODE 

Qf'' = litqf'''-iqf''?+e), q"^^ = x . (4.1) 

In particular, we will see that q^ falls into the right-hand well by a time of order y^e|liicT| 
after the bifurcation and so too does qt- The strategy is as follows: 

(1) Show that qt is of order y/e when t = y/e. 

(2) Show that {qt-, t) then leaves the space-time set K,{k) (see ()4.8p ). by a time of order 
yejlnaf. 

(3) Show that qt approaches the right-hand well and stays in a small neighbourhood 
of it up until any time ^2 > 0. 

(4) Show that by taking t2 large enough, qt stays in a neighbourhood of the right-hand 
well of order t^/2-7 for any < 7 < 1/2 and all t ^ ^2- 

Step One: We begin by describing how q^ behaves. 

Lemma 4.3. Let qf^*' he the solution of ^.1^ . Then we have, uniformly for all initial 
conditions x £ X , 

det^U/\t\for -T + e\lne\ ^t^ -^ 
\^/e for - ^/e i^ t i^ ^/£ 

and there is a constant C > such that for all x £ X , \qt\ < C for all —T ^ t ^ —T + 

e| lne|. 

Proof. Consider the equation 

tql{t)-q\{tf + e = ^. (4.2) 

By fixing t and differentiating the left-hand side with respect to q*^_^ we see that for i < 
it has no turning points and so admits a unique real-valued solution. Furthermore, we can 
check that qX{t) x e and {q\)' {t) = q\{t) / {^q\{tY — t) ^ e for negative t bounded away 
from zero. 

Suppose first that the initial condition satisfies x ^ q\{—T). Define zt = qf'^^ — qX{t). 
As long as zt ^ 0, we have zt ^ tzt/e so that 

^ gf * - qX{t) ^ (x - ql{-T)) e(*'-^')/2- . 

Let io = — T + e| lne|. If zt < for some —T < t < to, which means that qt ^ ^i then the 
analysis below for t ^ to can be applied from that time. Otherwise, the above inequality 
shows that qf^^ x e. 
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If X ^ 9+(— r) then we define zt = q\{t) — qf'^^. As (g!j_)'(t) > for negative t bounded 
away from zero, we have zt ^ for such t. In this case, there is ci > 0, independent of 
X & X, such that 

Zt i^cie + -{tzt + ?,qX{t)zi) 

for all —T ^ t ^ to- Furthermore, as long as zt ^ — t/6g^(t) (which is satisfied by z{—T) 
for all X ^ X by taking e sufficiently small) then ?>q*^{t)z1 ^ — tzt/2 and so 

1 

Zt < Cl£ + —tZt . 

This tells us 

^ q*+{t) - gf * ^ {q*{-T) - x) e(*'-^')/4^ + ciS [ e^*'"^')/^^ ds , 

J-T 
which shows that Zt ^ — t/6q*^{t) for all —T ^ t ^ to and we can use the above inequality 
to again see that qf^^ x e. 

We now analyse the behaviour for t ^ to- As qf^^ > and qf'^^ = 1 whenever qf'^^ = 0, 
it follows that qt ^ for all t ^ Iq. Therefore, for t ^ to we have 

qf^'^^itq^+e) 



and so 



^det 



Qt ^ Qto e 



* - - fc2e/\t\ for to ^ t ^ - Ve 



^ ^det^(.^-*g)/2. ^ / e(*^-^)/2-ds ^ "'^/l^' iu. .0 ^. ^ - V^ (4 3) 



|c2Ve for - Ve ^ * ^ \/^ 



for some constant C2 > 0. To obtain the lower bound, we use that for t ^ 0, tq — q^ '^ 2tq 
as long as q^ ^ |t|. By taking e sufficiently small, we have < qt^ ^ Vrol for all initial 
conditions x G A'. As long as ^ qf^^ ^ vN"' then 



(?f*^^(2tgf*+£) 



Jtn 



'to 
^det 



and 



By (I4.3P , we certainly have qt ^ Vrl for t ^ — ^/e and so the above inequality gives the 
corresponding lower bound for g^ up until this time. For —\/e ^ t ^ -y/e, we then have 
^gdct _ (gdct^s ^ _ (j^ £qj. gQj^g constant C > 0, so that g^*'* remains of order ^/£ in this 
interval. This completes the proof. 

D 

We now show that the deviation process yt := qt — qf'^^ satisfies |y(-v/e)| < he~^'^ for 
some h <^ e^'^, which will complete Step One. The process yt solves 

dyt = -Ht)yt + b{yt, t)] dt + ^ dWt , y{-T) = , (4.4) 

£ ye 

where a{t) = t — 3{qf^^)'^ and b{yt,t) = —Sqf^^yt — yf. For ah pairs {y,t) G B{h) for a 
choice of /i = ©(e^/'^) (see (gZl) and LemmaHUD, we have |6(y,t)| ^ My^. Solving (fO]) 
gives 

yt = ^ f e°(*'^)/^ dW^, + - /* e°(*'^)/^ 6(y„ .) d. =: y" + y] , (4.5) 
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where a{t, s) = J^ a(u) du. 

We now define the space-time set 13(h) mentioned above. If we write Var(yj') = a'^v{t), 
then we find that v{t) solves the ODE 

ev = 2a{t)v + 1 , v{-T) = . 

Let ^(t) be a particular solution of this ODE with nonzero initial condition, given by 

m = a-T) e2"(*'-^)/- +- f e2"(*'^)/- d. , a-T) = ^1 /^,| . (4.6) 

e j_rp z\a[—i )\ 

Then we define 

B{h) = {{y,t):-T^t^^,\y\<hy^} (4.7) 

and the stopping time tb(/i) = inf{t ^ — T : {yt,t) ^ 13(h)}. Before estimating rg(^), we 
must first understand how a(t) and ^(t) behave: 

Lemma 4.4. Let qf^^ solve ^4-1^ , define a(t) = t — 'i(qt^^)'^ and let ^(t) be given by |^.6p . 
Then, uniformly for x £ X, a(t) x t for —T ^ i ^ — y/e and \a(t)\ = 0(y/£) for \t\ ^ y/e. 
We also have, uniformly for x £ X, 

and \i(t)\ = 0(l/e) for all -T ^ t ^ ^. 

Proof. The assertions about a(t) follow from Lemma r4.3i We can use this to tell us how ^(t) 
behaves, for which it is helpful to consider the case a(t) = i as an example. Furthermore, 
since i solves e^ = 2a(t)i + 1, this tells us that |^(t)| = 0{l/e). D 

Having established the behaviour of all relevant quantities, we can now prove the fol- 
lowing proposition telling us that sample paths are likely to remain in B{h) for all times 

Proposition 4.5. There exists a constant C > such that for all e sufficiently small, all 
a < h <^ e'^''^ and all initial conditions x £ X , 

C { h"^ 1 

^{-^B(fe) < \/e} ^ ^exp<^-^(l-r(/i,e))^ , 

where r(h,e) = 0(^/£) + 0(he~^'^) uniformly for x £ X. 



-T,0| 



Remark 4.6. Choosing h = kay^\ lna\ with k > large enough guarantees that the right- 
hand side tends to zero as a I and that /ly^tv^ ^ V^> ^"^ which case we may take 
q(y/s) ^ -y/e uniformly for all x £ X . 

Proof. Recall the decomposition yt = Vt + Vt from (14. 5p . We have for all t < T^(h) ^ \f^-, 

\v\\ / 1 1 '■^ 



, ^ -^^ - / e"(*'")/^ \KVuM\ dtx 



Mh^ / \ 1 rt 

y/^ K-T ^ « ^ t '' ' V £ J-T 
ciMh^ 



^ ^!^ ( sup a^)) - I e"^*''')/^ du 



^ p3/4 
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for some constant ci > 0, where we obtain the final inequahty by bounding 

eJ-T " \c/^ioi\t\^^ 

and 



sup ^(u) ^ ' 



Tew V-T^«^t 7 " \C7e-i/4 for |t| ^ 7i 

Therefore, if |2/°|/\/^ < h{l - ciMhe'^/'^) for all -T ^ t ^ ^e then we must have 
T^ih) > V^- Letting H = /i(l — ciMhe~^''^), we obtain exactly as in the proof of Proposi- 
tion 4.3 from |BG02| that for sufficiently small e, 

for some C > 0. Note we cannot apply that proposition directly because our function 
a{t) behaves differently for |t| ^ ^/e than the corresponding function there. In particular, 
here a{t) < for t ^ e, whereas in [BG02J a{t) = t + 0{t^). However, in our case, 
\a{t, s)\ = 0{{t — s)^/e) for —^/e ^ s ^ t ^ -^e and we can check that this still allows us 
to suitably bound the term Pk appearing in equation (4.23) of their proof. 

D 

Step Two: We define for k > the space-time set 

/C(k) = {{q,t):t^V^,q^^{l- K)t} . (4.8) 



The boundary of /C(k) consists of the curves (ib-y/(l — K)t,t). For the present case, we 
only need to consider q ^ (for the slow pulling regime in Section |321 we will also consider 
q < 0). Let to ^ \/e and suppose that < q(to) < y^(1 — K)to. Then for t > to and as 



(4.9) 



long as 


0^ 


qt ^ ^(1 - ^ 


:)t, 


we have 


Qt ^ Qt: where 










dq^ 


= 


Ik, cr 

-Ktgf dt + —^ 


dWt, 


9"(to) 


= Q{to)/2 . 


Solving 


this SDE gives 




















Q^ = 


1 
2 


q{to)e^'"- 


-4)/2e 


a 


Jto 


e(*^- 


^')/2^ dWs . 



We now state two lemmas that are analogues of Lemmas 3.2.10 and 3.2.11 from |BG06| 
and are proved similarly. For the present section, we will only need to take to = \/£ and, 
by Step One, go = q{V^) ^ V^- For the slow pulling regime, other initial conditions will 
be considered. Let t;(;(k) = inf{t ^ to : {qt,t) ^ ^('«)} and r^ = inf{t ^ to : (/f ^ 0}. 

Lemma 4.7. Assume that qt starts at time to ^ ^/e in go > 0, where (goi^o) £ JC{k). 
Then there is C > 0, independent of to, qo and k, such that for all t ^ to + s/to, 

a 
Lemma 4.8. Let gf start at time to ^ \/£ in go/2 > 0. Then there exist constants 
Ci, C2 > 0, independent of to, go and k, such that for all t ^ to, the probability of 
reaching zero before time t satisfies the bound 

P*«'«»K <t}^ -^exp|-^^l . 
QoVto I cr^ J 
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This second lemma shows that when (/oV^ ^ '^i ^^^ hnear process gf is unhkely to 
return to zero for any time t ^ to, while the first lemma shows that as t increases the 
probability of qt remaining in /C(k) decreases. Therefore, we have qt ^ (?f until time 
"OcCk) ^iid so qt must exit /C(k) through the curve ^J{1 — K,)t. Indeed, there are constants 
Ci, C2 > such that for initial time to = ^/^ and any initial position go >< V^-, we have 
for any t ^ 2^/e that 

IP^'-{r^W < t, r« > t} ^ 1 - ^V^e--^/^^ " I^^^P {-^} ' ^^-^O) 

In the present fast pulling regime, e ^ a'^'^ and so the third term on the right-hand side 
tends to zero as u ^ 0. Picking t = y^2te|lria| , we see the second term also tends to zero 
as long as A; > 1/k. By a time of order ^Je\ lncr|, all paths will have left /C(k) through its 
upper boundary. 

Step Three: Firstly, we will see how deterministic solutions behave when started from the 
boundary of /C(k). For this, we let q\{t) be the same solution of (j4.2p that we considered in 
the proof of Lemma 14.31 i.e. the unique real-valued solution existing for all times t ^ —T. 
For t > and e sufficiently small, we have \/t ^ Q+{t) ^ \/t -|- e/t, which can easily be 
seen by the intermediate value theorem. For the sake of brevity, we shall write r to mean 
tjc{k) ■ The following proposition tells us how q^ ^ '^ behaves, where q^ ^ '^ solves 

'Zf *'^ = -^{tqf^''^ - (^^'0' + e) , 9^'^ = a/O^^ . (4.11) 

Proposition 4.9. Assume that k G (1/2,2/3) and let ij = 2 — 3k, > 0. There is a constant 
C > such that the solution, qt^ ''^ , of |-^.ii| ) satisfies 

^ qlit) - gf '^ ^c[-^ + {q*(T) - q'^^^'^r)) ^-^^^'-^')l^- 

for all t ^ T and e sufficiently small. 

Remark 4.10. The condition k > 1/2 guarantees that paths do not re-enter IC{k) after 
leaving, while k < 2/3 ensures that the potential is convex outside of IC{k). 

Proof. The inequality qt'^ '^ ^ Q+i^) follows since g'^''*''^(r) < (7+(t) and 

The proof of the other inequality follows along the same lines as that given in [BG021 
Proposition 4.11]. Note, however, that unlike there we only need to take e sufficiently small 
and not t. This is because in our case the value of Oq, which is defined in equation (4.99) of 
|BG02j . is given by -2(1 + 0^(1)), rather than -2(1 + 0^(1)). Similarly, M* = 3(1 + 0^(1)). 
As q\{t) ^ Vt + e/t < (3/2) \/t for t ^ \/e and e small, we can use (14.121) to show 
(g;)'(t) ^ (3/4)t-i/2, giving K* = 3/4, where K* is also defined in (4.99). In [BG02], 
K* = 1/2, but the proof just requires that K* < 1. D 

Now that we understand how q^ '^ behaves, the final step is to show that qt, starting 
at the same point, stays close. Having shown that the analogue of Proposition 4.11 from 



16 MICHAEL ALLMAN, VOLKER BETZ AND MARTIN HAIRER 

|BG02j holds, the proofs of the subsequent bounds there can easily be extended to our 
case and we now show what these are. Let 

where a^{T) = t—3{q^^ '^)^ is the linearisation of the drift term around q^^ '^ and a'^{t, s) = 
f a{u) du. As is shown in Lemma 4.12 from |BG02] . it follows from Proposition 14.91 above 
that |a^(r)| x i so that ^^(t) x 1/t. 
Now we write 

A^{h) = {{q,t):t^T,\q- qf^'n ^ h./W)} (4-13) 

and let T4t(/j) = inf{t ^ r : {qt,t) ^ A^{h)}. The following bound on t_4t(/j) follows by the 
analogue of Theorem 2.12 in [B G02] . It tells us that for k £ (1/2,2/3) and any t2 > 0, 
there exist constants C, ho > such that for h < /iqt and e sufficiently small, 

p.,Va^;^{r^.(,)<t2}^^exp|-i^[l-0(.)-o(^)j| . (4.14) 



The right-hand side becomes small by choosing h = kay^\ lna\ for k large enough, for 
which we note that r ^ ^/e by definition so that h <^ r. 

Step Four: Let us suppose that t2 ^ 1. By Step Three, we may write q{t2) = \/t2 + 
0(a^\ ln(T|) + C'(e) independently oiT^^i^^-^. For a given g(i2), let g^*^* be the corresponding 
deterministic solution starting at the same point. We can again obtain a similar bound 
as in Proposition 14.91 but for simplicity let us just say that (9/10)\/t ^ (7^°' ^ (ll/10)\/i 
for all t ^ 1. Letting yt = qt — (if^^, we define t(7) = inf{i ^ ^2 : \yt\ > t^'"^'^} for 
< 7 < 1/2. We again decompose yt into a linear part, y^, and nonlinear part, yj, as 
in (j4.5p . Then a{t) = t — 3(gf°*)^ x —t uniformly for all t ^ 1 and the function b{y,t) 
containing the nonlinear terms now satisfies \b{yt,t)\ < My/tyt for all t < t^j) and some 
constant M > independent of t2- We will show that P{r(7) < 00} — t- 0. For t ^ t(7), 
we have 



y}\ 


^ 


i/'e 


a{t,. 


)/e 


\h{yu 


u)\ du 




< 


e 


-27 


f 

Jt2 


e-(M 


)/" du 



where the final inequality holds uniformly in t and the constant C > is independent of 
t2. Therefore, if |y?| < H{t) for all t ^ t2, where H{t) = t^/2-7(i _ Ct^^), then we must 
have r(7) = 00. Note that 1 — Ct~'^ > and H{t) > for all t ^ t2 by taking ^2 large 
enough. We have 

-p||>iUf;4 -p ||>4' 

t>t2^[t) J ^ [s, sgt^Sj + i -H(i) J 

where ^2 = so < -si < . . . is chosen by —a{sj+i,Sj) = e^. Note that, uniformly in j, we 
have s'j_^_l — Sq x — q(sj+i,so) = — a(sj-|_i, Sj) — ... — a{si,so) = (j + l)e^, which shows 
that Sj — )• cxD as j — )• 00. Call the summand on the right-hand side above Pj. As H(t) is 
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increasing, we can further bound Pj by replacing H{t) with H{sj). We can also use for 



Sj ^ t ^ Sj+i the inequality 

t 



Wt 



Ol 



^ / e"(*'^)/^ dW. 

V^ Jso 

This gives for all j ^ 0, 

ft 



Pj^ 



sup 

So ^ t ^ Sj + 1 



M^)/^ dW. 



^ / e-°(^)/" dWs 

V^ Jso 
,}11q-'^{so)/^H{sj) 



^ 2 exp '' — 



2cj2 P+^ f,2a{sj+^,s)/e ^ 

J SI) 



^ 2 exp 



ciSjH{sj)'^ 
2^2 



where the constant ci > in the final inequality is independent of j. Note that the second 
inequality comes from Lemma B.1.3 in the Appendix of |BG06j and the final inequality 



uses that a(sj+i,s) 



j^i — s ) uniformly for all sq ^ ^ ^ ^j+i ^^'^ ^11 J- Summing 



over j ^ 1 and using that Sj — Sj-i ^ Ce /sj uniformly in j, we have 



E^. = E^. 



i=i 



i=i 
C 

7^ 



So Sj — l 



^ 2 A^Pj^n^J ^j-'i-j 

3=1 
/•oo 

^C sP{s)ds, 
Ju 



where 



P{s) = 2 exp 



'i2 

cisH{sf 
2^2— 



^ 2 exp 



C2S 



2(1-7) 



cr^ 



and C2 > is a constant. We have used that sP{s) is decreasing when bounding the series 
by the integral above. Then 



/ sP{s)ds^Ca exp{— C2t 

Jt2 



2(1-7) 



for some constant C > depending on ^2 and 7, so that 



> U ^Pq + 



and the right-hand side tends to zero as o" J, 0. 



Ca' 



/-'} 



2(1-7) 



exp{-C2t2^^-^7<^'} 



4.2. Slow pulling. We now consider the slow pulling regime from Proposition 14. 2i In 
this case, the noise dominates the dynamics and cancels out the asymmetry caused by 
pulling. The process qt should, therefore, behave similarly to qt, where 



1 



a 



dqt = - {tqt - qf) dt + ^ dWt , q{-T) = . 

£ \J E 



(4.15) 



As we have chosen q{—T) = 0, the law of q is entirely symmetric about zero. The strategy 
is as follows: 
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(1) Recall from |BG021 IBG06J that qt stays close to the origin with high probability. 
At time t = y^, its typical spreading is of order ae~^^^^/\ lna\. 

(2) Show that paths of qt stay close to those of qt until qt leaves the diffusion-dominated 
strip S{h) defined below. 

(3) Show that qt then exits the slightly larger strip, /C(k), without returning to the 
origin. 

(4) Show that qt then finally falls into the potential well on the same side as it left 
IC{k) and remains there. 

Step One: This step is the same as Step One from the previous section, except now we 
are analysing the behaviour of qt- Unlike in the previous section, we can use directly the 
results of |BG021 IBG06) . which we now summarise. We again define the function ^(i) as 
in the last section and now a{t) := t — 3(gj )^, where 

fit' = \mf'' - {qt'f) , e^ = . 

Clearly, q^""^ = and so now a{t) = t. Again, ^(t) x l/(|i| A ^) for -T ^t ^ ^. We 
define the space-time domain 



B{h) = {{q,t) -.-T^t^^e, \q\ < /ly^} 

and the stopping time Tg(;i) = inf{t ^ — T : {qt,t) ^ t3{h)}. Applying Theorem 2.10 from 
[BG02J . we see that there exist constants C, /iq > such that for e sufficiently small and 
h ^ hoy/e, 



Choosing h = kay^\ \na\ for k large enough, the right-hand side tends to zero. At time 
^/e, we may take |(7(\/e)| = 0{ae~^^'^^/\]na\). 

Step Two: In the fast pulling section, we saw that at time ^/e, q{\/s) x \/£, from which 
we could then show its subsequent exit from /C(k). Now we are considering the exit of qt 
from /C(k) and with other values of q{\/s) as found in Step One. Before looking at the exit 
of qt from /C(k), we must first analyse its exit from a smaller strip. We define for t ^ ^/e 
the diffusion-dominated strip 



S{h) = \{q,t):t'^^,\q\ < A| 



and the stopping time t^^^j) = inf{t ^ y/e : {qt,t) ^ <S{h)}. See Figure 3 in fBG 02| and 
Figure 3.12 in [BG06J for an illustration of 5(/i) and /C(k) (note that IC{k) is denoted 'D{k) 
in |BG02| ). Let h* := /iq cry^\ lna\, where /iq > is a constant sufficiently large so that 
{q{^/e), \/e) G S{h*). Applying Proposition 4.7 from |BG02] with the choices h = h* and 
^ = 2, we see that there exists C > such that for all a sufficiently small and all initial 
conditions {qo,\/£) G <S{h*), 

^^''"'{rsm^t}^ (^) exp|-^^^^[l-0(l/ln(r/a))]} , 

as long as o"| Incp^ = 0{y/e), which we already assume in the slow pulling regime. We 
can check that by taking t = ^J2kehi{h* /a) with A; > sufficiently large, the right-hand 
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side tends to zero. For such a choice of k, we define t* = y^2keln{h* /a) and henceforth 
assume t^j-^j.) ^ t*. 

The important point here is that by symmetry, qt exits S{h*) through either boundary 
with equal probabihty. Now that we understand the behaviour of qt up until its exit from 
S{h*), we turn to qt- The following lemma shows that qt is close to qt at time r^^/j.). 

Lemma 4.11. If s{a) < a^/^| lncr|"^^/^ and T5(/i*) ^ t* , then 

/ ,, /e3/4|i^^|13/8\\ 
Qrs,,*) = qrs.nn I 1 + O I 1 1 • 

The proof of this lemma is based on the following simple comparison of qt and qt- 

Lemma 4.12. Let qt solve \2.4^ with initial condition q{—T) = x €z X and let qt solve 
l^.iJp . We have, almost surely, for all t ^ — T, 

ft 
Qt 



if X ^0 and 

ifx^O. 
Proof. We have 



J-T 



_1 /■* e(*'-^')/2-(?i'ds + ^ /* e(*'-^')/2- dVF, . (4.18) 
e J-T v^ J-T 

By the comparison principle, qt ^ (Jt almost surely, where qt solves 

dqt = -{tqt-qf)dt+^dWt, q{-T) = x . 

£ ye 

Using this lower bound for q in the second integral on the right-hand side of (j4.18p gives 

qt^qt + f e(*'-')/2^ ds . 

When j; ^ 0, we have qt ^ qt almost surely, which gives the upper bound in (j4.16p . For 
the lower bound, we have 

g, = xe(*'-^')/2- -- f e(*'-^')/2-g3ds + ^ r e(*'-^')/2- dl^, 
e J-T v^ J-T 



£ J-T V^ J-T 



The case x ^ is easier and does not involve qt- It follows along similar lines. D 



20 MICHAEL ALLMAN, VOLKER BETZ AND MARTIN HAIRER 

Proof of Lemma \4-ll ■ For all ^/e ^ t ^ t* , we have 



r e^^'-^'y^'ds^CiV^f^Y ^C2V^\lnaf (4.19) 



and for all x G A", 

4 



Of these two estimates, (j4.19p gives the larger upper bound. Next observe that 



Therefore, 



V^sm < ^_ e^/^ 



Using (j4.19p and the above inequality together with Lemma 14.121 gives the result. D 

Step Three: Now we analyse the behaviour of qt, rather than qt, for t > T^{h*)- ^.s we 
saw in the fast pulling case, if qt starts at time to ^ V^ at q^ > with ((^Oi^o) ^ IC{k) 
then qt ^ gf as long as < qt < y/{l — K)t, where qt was defined in ()4.9p . Unlike in the 
previous section, we now have to also consider negative initial conditions. We would like 
a bound of the form qt ^ qf in such cases, but now the bias of qt in the positive direction 
makes this more difficult. In order to obtain a corresponding comparison with q^, we need 
an additional assumption on to a-nd qo as set out in the following lemma. 

Lemma 4.13. Suppose that at time to ^ \/e, Qt starts at qo < 0, where {qo,to) G JC{k) and 
\qo\ ^ e/to- Then we have qt ^ qt o,s long as —^/{l — n)t ^ qt ^ and e is sufficiently 
small. 



Proof. For — y^(l — K)tQ < go < 0, it is certainly true by comparison of the drift and initial 
conditions that qt is bounded above by solutions of 

dzf = -{ntz^ + e) dt + ^ dWt , ^ = go , 

e ye 

as long as — y^(l — K)t ^ qt ^ 0- The result follows since 

Z- = goe'^(*'-*o)/2^ ^ /■* e«(t2-s2)/2e^g^ ^ /■* g.(t2_s2)/2e ^^^ 
J to V £ J to 

^ (go + Ce/to) e''(*'-*o)/2^ + ^ r e'^(*'-^')/2. ^^^ 

n 

This lemma shows that, under suitable conditions, we may compare qt and gf for both 
postive and negative initial conditions go. For go < we get analogous bounds to those in 
Lemmas 14.71 and 14.81 In the previous section, we applied those lemmas with to = y/e and 
qo ^ V^- We now apply these lemmas with to = Tg{h*) and |go| x h* /^/t^(J^ (see Lemma 
I4.1ip . Note that if \/e ^ Ts{h*) ^ i* then e/to ^ |go| and so the conditions of Lemma 
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14.131 are satisfied. We obtain a bound similar to (j4.10p and again see that IC{k) is left by 



a 



time of order y/e\ lna\. 



Step Four: When qt exits /C(k) on the positive side, this part is exactly the same as Step 
Three from the previous section. The other case when qt exits /C(k) on the negative side 
is similar. Firstly, we introduce gl(i), another real- valued solution of (|4.2|) existing for 
t ^ ^/e and satisfying the bounds — \/f ^ ql{t) ^ — \/t + e/t and (gl)'(t) < for all such 
t and £ small. Note there is also a third real- valued solution of (j4.2p between q!l and g!^, 
which is an unstable equilibrium branch. By taking e small, we have g!L(t) < — y(l--re)i 
for all t ^ ^/e. Again writing r = Tj(^(k)j ^^ need to check that the deterministic solution, 
q^'^ '^, of (|4.1ip . with initial condition ^'^'^^''^(t) = —y^T^^l^, satisfies a bound as in 
Proposition 14.91 of the form 



^ q^'^ - q*_{t) ^ C ^-^ + (g-^'-(r) - g*(r)) e" 

for all t ^ T. The main thing to ensure is that q^^ '^, which has a bias in the positive 
direction, does not re-enter the set IC{n) after having left. To see that this is indeed the 
case, first note that the derivative with respect to t of the boundary curve, —^/{l — K)t, is 
given by — 2t~^'^\/T^-^- The derivative of q^'^ '^ when on the boundary of IC{k) is given 
by jt^''^{—K + e)\/l — K. Using that t ^ -y/e, we see that the inequality 



-t^/^(-K + £)Vi -K< --r^/Vi-K 

holds when k > 1/2 + e, which is true for all k > 1/2 by taking e sufficiently small. Having 
established that qZ. (t) ^ 9^ ^ '^ ^ — -y/(l — K)t for all t ^ r and e sufficiently small, the rest 
of the proof follows like Proposition 4.11 from [BG02J . The subsequent estimate ()4.14p 
above showing the concentration of qt in the set ^^{h) then follows, where ^^{h) was 
defined in ()4.13p . Finally, we can show as in Step Four from the fast pulling section that 
qt stays in a neighbourhood of — \/t for all t ^ t2. 

5. The Full Solution 

We now consider the full equation ()2.2p and show that for sufficiently small mass (large 
/3), it behaves like the overdamped solution of the previous section. The general strategy 
is as in Section U namely we first show that if we start the system at the origin at time 
s <C — 1, then the solution at time —1 belongs to a suitable set. This is done in the 
following two propositions. We then provide a result that is uniform over all solutions 
starting from the set in question. The first step is achieved by the following statement: 

1_fl 1_fli ~ l+3,a 1+3/3 

Proposition 5.1. Let X = [—e P,e '^], V = [—e 2 ,e 2 ], letT '^ 1 he a constant, 
and let qt solve \2. 2\) with f3 > 2. Then, we have 

lim liminfP*(g_2r G ^ , P-2T G V) = 1 • 

CT,e— s>0 s— >-oo ^ ' 

Proof. We fix some arbitrary starting time s < —2T and we consider the solution to (j2.2p 
with initial condition qs = Ps = 0. We define the function ^{p, q, t) by 

e^ i 1 1 

^{p, q, t) = yp^ - ^9^ + -^q^ -q + ^pq- 
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Applying Ito's formula to "^{pt,qt,t), we obtain 

d^{pt,qut) ^ (-e-^^{pt, qt, t) + \<^^e-^~^ " ^'z' - ^% j dt + dM{t) , 
where M is some continuous martingale. Using e^~^q'^ + <? ^ — (l/4)e^~^, we see that 

^E*(pi, qt, t)^ - e-^E^{pt, qt, t) + ^^'e-'-" + ^e'-^^^ . 
It follows immediately that E^{pt,qt,t) ^ o-^/(2e) + {l/8)e^^'^ . Since for t ^ - 2 we have 

it follows that 

Eq^ ^ Se^ + £^-2/3 + 4^2 ^ ^2-/3 ^ 

Ep2 ^ 4^1-/3 + (i/2)ei-3/3 + 2^2^-1-/3 ^ (i/2)£l-2/3 . 
The stated result then follows at once from Chebychev's inequality. D 

Now we use the previous proposition to restart the process at time — 2T. We denote by 
T^x,v ^YiQ law of the solution of (j2.2p starting at time —2T with q-2T = x, p-2T = v. 

Proposition 5.2. Let X = [—1, 1], V = [—£~^,£~^], let T ^ 1 be a constant, and let qt 
solve /i2. 2\) with (3 > 2. Then, we have 

lim inf P^'*(g_T G Af , p_T G V) = 1 . 

1 e/^ £2/3 g/3-1 

^{P, Q, t) = p^ + yP? + ^/ + ^/ • 
Applying Ito's formula to ^{pt,qt,t), we obtain 

1 / e^^'^ tie a'^\ 

d^ipt, qt, t) = - ( —pi + -qi - -qf + -qt + eHptqt + e^+V + y ) ^^ + d^^(*) ' 

where M is some continuous martingale. Using that t G [— 2r, —T] and /3 > 2, we have 
for sufficiently small e, 

d^(pi,gi,t) ^ i (-^*(pt,%,t) + 4e'^+' + I + y) dt + dM(t) . 

It follows that E^(pi, gt, t) ^ e-(*+2^)/2e E*(p_2T, g-2T, -2T) + a^ + Se'^+i + e. We can 
then use the bounds '^{p,q,t) ^ e'^^p'^/4 and '^{p,q,t) ^ g^/8 to obtain the result. 

D 

We would like to use a singular perturbation approach to show that for t ^ — T and 
suitably large /3, sample paths of qt can be approximated by those of an overdamped 
equation starting at —T. For this, we need to first consider a process that is similar to qt 
but has better regularity. To this end, we introduce the processes Q and P that solve 

dQt = Ptdt, Q_T = P-T = 0, 

e^dPt = -Ptdt + ^dWt. ^^■'^^ 

Ve 



Proof. Let 
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We find that 



and 





Pt 


= ae 


1/2-/3 


L 


e 


-{t~s)e- 


-/3 


dWs 




Wt- 


e^Pt 


= ae~ 


1/2- 


-/9 


/>- 


-{* 


-s)e-^ 



Q, = --W^j - e^p^ = ae-'/'-" / e"^*"^^^ W{s) ds 

For 5 > we now define two events, Ei and £'2, by 

El = {\Qt\ > ae-^/'^-^ for some t G [-r,i2]} , 
E2 = {\Pt\ > fje-^/2-/3/2-5 fQ^ gQj^g ^ ^ [-T,t2]} . 

Lemma 5.3. There exists C > such that for all 6 > and all ^2 > 0, 

¥{EiUE2)^Ct2{e-^-'''''^^' +e--"/2) 
holds for a, e > sufficiently small. 

Proof Since Qt = ae-^/^-/? j*^ e-^*-*)^"'' W{s) ds, it follows that 
F{Ei) ^ P{|W^i| > e-^ for some t £ [-T,t2]} 
€ Cto e" 



_^-S/{t2+T) 



where C > is independent of 5 and ^2 • We also find that there exists C > independent 
of 6 and ^2 such that for e sufficiently small, 

¥{E2) ^ Ct2e-^^-^ exp |--e-2'5 

The result follows by combining these two estimates. D 

Now let yt = Qt — Qt- It solves, almost surely, the second-order ODE 

e^y = -y + ^{t{yt + Qt) - {yt + Qtf + e) , y{-T) = x , y{-T) = V . 

The following proposition shows that for almost all paths in {Ei U E2y, y may be ap- 
proximated by the solution of a first-order ODE. Note that the condition on /? is a little 
stronger than necessary, but is required later on in this section. 

Proposition 5.4. For all ^2 > 1; there exists C = C{t2) > such that for all f3 > 2, all 
<5 < (3/2-1, all 0-2/(^+2^) < e < 1 and almost all paths in (Ei U E2y, 



1 „ 

y - -{t{yt + Qt) - {yt + Qt)^ + e) 



^Cmax{e^-2-^ae-3/2+/5/2-25| 



for all t G [— T + 2e^ ,^2] O'^d a sufficiently small. 

Proof. If we write z = y, then almost surely the pair (y, z), which are differentiable, solve 

y = z , 



e^z = -z + -g{t,yt + Q 



1 



e 

where g{t, yt + Qt) = t{yt + Qt) - {yt + Qt)^ + £■ For t > 0, we have g{t, y/t) « so that 
we do not expect yt + Qt, or indeed yt, to be much larger than y/t. Therefore, we let 



24 MICHAEL ALLMAN, VOLKER BETZ AND MARTIN HAIRER 

r = inf{t ^ — T : \yt\ > 2-y/t^}. On {Ei U £'2)'^, there is C > depending on ^2 such that 
for all —T ^ t ^ T At2, \g{t, yt + Qt)\ < C. We solve the equation for z to give 



zt = v e-'-'+n^-' + e-(i+W / e-^*-^^^"^ g{s, y, + Q,) ds (5.2) 

almost surely, from which we deduce that l-^tl ^ g-P Q-\t+T)e _j_ (jj^ £qj. all t ^ r A ^2- 
This immediately shows that for —T ^ t ^ {—T + 2e^~ ) A r and sufficiently small e, 

|yt - xK 1 + Ce'^-'^-^ (5.3) 

and SOT > -T + 2el^-^. 

For — T + e^~ ^ t ^ r A ^2) we have l-^il < C/e. Furthermore, for such t we find 

d 



-^g{t,yt + Qt] 



^Cmax{E-\ae-'/'-^/'-'}. 



Now we apply the Laplace method to the integral in ()5.2p . For t ^ — T+2e'^ , decompose 
the integral as 

rt rt-e^'^ ft 

+ / 
Then, by the boundedness of 5, we have 

rt-eP-^ 



j ' e-(*-^)^"%(s,y, + Q,)ds 



<Ce-^"* 



For the remaining integral, we use a Taylor expansion of g to give 

g{s, Vs + Qs) ^ g{t, yt + Qt) + C{t - s) max{e~\ae-^/^-^/^-^} . 
Then 

ft 



[ e-^*-^)^"" g{s, Vs + Qs) ds ^ e^g{t, yt + Qt)+ 

+ CePe-'-' +Cmax{e2/3-i-5,^e-i/2+3/J/2-25| ^ 
which tells us that for —T + 2e"~ ^ t ^ t At2 and u sufficiently small, 

^* ^ lg{t,yt + Qt) + Cmax{e/^-2-5^^^-3/2+^/2-25| (54^ 

In a similar way, we can also show that 

zt ^ lg{t,yt + Qt) - Cmax{e/^-2-^ae-3/2+/^/2-25} . (5.5) 

We will now show that the assumption r ^ ^2 leads to a contradiction. For this, we will 
show that if yr = +2y/t2, then the right-hand side of (j5.4p is strictly negative, whereas we 
should have z,- ^ by continuity. The case yr = —2\/t2 is similar. First, we note that if 
y^ = +2y/t2, then there are constants Ci, C2 > depending on ^2 such that 

g{r,yr + Qr)^ -Ci + C2ae-'/^-^ 

and for a sufficiently small, the right-hand side is strictly negative and bounded away 
from zero. Then the conditions on e, (3 and S guarantee that the right-hand side of ()5.4|] 
is strictly negative. This means that we must have r > ^2 and so (15. 4p and ()5.5p hold for 
all — T + 2e"~ ^ t ^ t2, from which the result follows. D 

Now we will use Proposition 15.41 to tell us something about the SDE (|2.2p . 
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Proposition 5.5. For all t2 > 1, there exists C = C{t2) > such that for all j3 > 2, all 
< 5 < /3/2 - 1, all 0-2/(1+25) ^e ^l and almost all paths in (Ei U £2)", 

qt ^qt + e^Pt ^ qt 
for all t G [— T + 2e^~ ,^2] cL^d a sufficiently small, where 

dqf = -{tqt - {qtf + e{l±r{a)))dt + ^dWt , q^i-T + 2e^-') = x ±3 , (5.6) 

£ ye 



with Wt the same Brownian motion appearing in \2. 2\/ and 

r{<y) = Cmax{e''-2-5^^^-3/2+/3/2-25| 

Proof of Theorem \5.5[ We will show the upper bound. The lower bound is similar. Letting 
to = —T + 2e"~ , we know by Proposition 15.41 that there exists C > and a process y^ 
solving 

yt = lait, Vt + Qt) + C maxje^-^-^ ^e-^l^^^l'-^^} , y+ = x + 2 , 

such that yt ^ yf for all t G [^0)^2], where yt = qt — Qt- Note the initial condition for y^ 
is chosen using (j5.3p . Therefore, qt ^ yt ~^ Qt- In the same way as in Proposition 15. 4^ we 
can show |y^| < 4-^/t^ for all to ^ ^ ^ ^2 by considering the sign of y^ ■ 

Now define r]t := yf + Qt + £^Pt and note that for all to ^ ^ ^ ^2, |^t| < C' for some 
constant C > depending on t2- It solves 

dr?t = -[t{r^^-ePp^)-[r^^-ePPtf+e+Cuis.^{eP-^-\ae-^'^+^/^'^^})dt+^dWt - (5.7) 
£ yS 

with initial position r/tg ^ x + 3. 

Denote the drift term above by f{t,r]t,Pt,£)- We will now show that / is bounded in 
such a way that allows us to use a comparison principle. As we are working on (Ei L)E2y, 
we know that \£^Pt\ ^ (Te"^/2+/3/2-5 ^^^ ^^l t £ [to,t2] by definition. Note also that by the 
conditions on 5 and e, max{e^^-'^^'',cre~^'2+/3/2-2<5| ^ ^^ -^^ then have 

fit, rit, Pt, e) ^ -^{trit - vl + e{l + r{a))) , 

where r{a) = Cmax{e'^~2-'5^(j£-3/2+/3/2-2<5| ^^^ gome constant C > depending on t2- 
Let qf' be the solution of 

dqt = -{tqt - iltf + ^(1 + ria))) dt + ^ d^ , q+=x + 3. 

£ \J £ 

By Lemma [A. II below, r/j ^ q^ for all t G [to)i2] and almost all paths in {E\ U E2)^. 
Therefore, qt ^ r]t — £^Pt ^ qt ~ ^^ Pt, which gives the upper bound. D 

This proposition allows us to complete the proof of Theorem 12.31 

Proof of Theorem \2.3[ Firstly, let us consider the fast pulling case. That is, cr^' ^ | In cj| ^/'^ <c 
e <C 1. In this case, when /3 > 2 the term 1 — r(o") appearing in (15. 6p is strictly positive and 
bounded away from zero for all a sufficiently small. Then the analysis of Section 14.11 can 
be applied to q^ (the slightly different drift term does not matter). By Lemma [5.3t we can 
assume {EiUE2y to hold, in which case \£^Pt\ ^ o-e~i/2+/5/2-<5 for ah t G [-r + 2e^-'^, t2]. 
Therefore, we find that q^ — £^Pt > 7\/t for all ciye|lna| ^ t ^ t2, where ci, 7 > are 
suitably chosen constants. 
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For the slow pulling case, let < (5 < (3/2 - 1 and cj2/(i+2'5) < e < (j^/^\\n(j\-^^l^ . 
Then again 1 ibr(o") is positive and bounded away from zero for a small and the analysis in 
Section W?2\ applies to qf and q'^ . Note that, in the limit, q^ and qf must "go the same way" 
by comparison of their drift terms and initial positions and we know by Proposition 14.21 
that the probabilities are 1/2 in either direction. As above, the term e^Pt does not change 
anything. Therefore, qt behaves in the same way as q^ and qf . D 

Appendix A. A comparison principle for SDE 

We present a lemma that is a slightly modified version of Theorems 5.1 and 5.2 appearing 
in |And72j . Our proof follows those given there. 

Let (Jl, J^, P) be a probability space on which is defined a one-dimensional Brownian 
motion W adapted to a filtration {J-t)t ^ o- For t ^ 0, let X[t) and Y{t) be two real-valued 
processes evolving according to 



X{t)=XQ+ I a{s,X{s),Z{s))<ls + CWt, 
Jo 

Yit) =xo+ [ c(s, Y{s)) ds + CWt , 



(A.l) 



where C is a constant, a : [0,oo) x M^ — )• M, c : [0,oo) x M — )• M are continuous functions, 
Z{t) is an T'i-adapted process with continuous sample paths almost surely and Z{0) = zq, 
which is J-Q-adapted. Note that X{t) and Y{t) have the same initial condition, which may 
be constant or random as long as it is J-Q-adapted. 

Lemma A.l. Suppose there are constants Ci, C2 > such that whenever \x\ < Ci and 
\z\ < C2, a{t,x,z) < c{t,x) for all t ^ 0. If, almost surely, \X{t)\ < Ci and \Z{t)\ < C2 
for all t ^ then 

¥{Y{t) ^ X{t) for allt^id} = 1 . 

Proof. Define r = inf{t > : Y{t) - X{t) < 0} and set r = +00 if Y{t) ^ X{t) for all 
t ^ 0. This is a stopping time because if t > then 

{r>t}= fl {Y{r)-X{r)^0}eTt. 

rG[0,i]nQ 

Put D = {t < +00} and assume that ¥{D) > 0. Then we can define a probability measure 
Q(-) =F{-\D) on J". Let 

X+{t)=X{t + T) , 

Y+it) = Y{t + T) . 

By continuity, Q{X^(0) = y^(0)} = 1. For any t ^ 0, we can therefore write (almost 
surely with respect to Q) 

Y-^{t) - X+{t) = Yit + r) - X{t + r) 

rt+T 



Y{t) - X(r) + / c{s, Y{s)) - a{s, X{s),Z{s)) ds 

t 

c{s + T, y+(s)) - a{s + T, X+{s), Z+{s)) ds . 
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The right-hand side is continuously difFerentiable in t and so 

y+(t)-x+(t) 



Q < hm ■ 



t 



c(r,X+(0))-a(r,X+(0),Z+(0)) 



1 . 



Therefore, Q{y^(t) > X'^i^t) for ah sufficiently small t > 0} = 1. But due to continuity 
of X and Y and the definition of r, this probability should be zero. This contradiction 
arises from the assumption that ¥{D) > 0. Therefore, P{r = +00} = 1. 

D 
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